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Soil tillage practices have a profound influence on the physical properties of soil and the greenhouse gas 
(GHG) balance. However there have been very few integrated studies on the emission of carbon dioxide 
(C0 2 ), methane (CH 4 ) and nitrous oxide (N 2 0) and soil biophysical and chemical characteristics under 
different soil management systems. We recorded a significantly higher net global warming potential under 
conventional tillage systems (26-31% higher than zero tillage systems). Crucially the 3-D soil pore network, 
imaged using X-ray Computed Tomography, modified by tillage played a significant role in the flux of C0 2 
and CH 4 . In contrast, N 2 0 flux was determined mainly by microbial biomass carbon and soil moisture 
content. Our work indicates that zero tillage could play a significant role in minimising emissions of GHGs 
from soils and contribute to efforts to mitigate against climate change. 

Globally, agriculture accounts for 10-12% of total anthropogenic emissions of greenhouse gases (GHGs), 
estimated to be 5.1-6.1 Gt C0 2 -eq yr~' in 2005 1 . Conservation tillage is one among many different 
mitigation options suggested to reduce GHG emissions from agriculture. Conservation tillage practices 
such as reduced/minimum/zero tillage, direct drilling and strip cropping are also widely recommended to protect 
soil against erosion and degradation of structure 2 , create greater aggregate stability 3,4 , increase soil organic matter 
content, enhance sequestration of carbon 5,6 , mitigate GHG emissions 7 and improve biological activity 8 . Derpsch 9 
estimated that approximately 45 million hectares was managed by conservation tillage worldwide in 2001 and this 
figure had more than doubled by 2007. 

Minimum tillage practices have been reported to reduce GHG emissions through decreased use of fossil fuels in 
field preparation and by increasing carbon sequestration in soil 10 . For example, Hermle et al. 13 observed net 
carbon sequestration to a depth of 50 cm after 20 years of no tillage. However, reduced tillage can lead to a 
stratification of soil organic carbon at the surface 11 in contrast to the more uniform distribution of carbon in 
conventionally tilled soils 12 . The crop residues accumulated on the soil surface under reduced tilled conditions 
may result in carbon being lost to the atmosphere upon decomposition 10 . Furthermore, climate change mitigation 
benefits such as reduced C0 2 emissions, by virtue of increased sequestration of carbon and increased CH 4 uptake 
under reduced tillage, could be offset by increased emissions of N 2 0, a greenhouse gas with higher warming 
potential than both C0 2 and CH 4 13 " 15 . Increased N 2 0 emissions have been linked to increased denitrification 
under reduced tillage due to the formation of micro -aggregates within macro-aggregates that create anaerobic 
micro sites 13 with increased microbial activity leading to greater competition for oxygen 16 . 

Reduction of tillage can also create increased soil densification and a subsequent decrease in the volume of 
macropores 17 leading to reduction in gaseous exchange. Soil aggregation and the resultant geometry of the pore 
structure are vitally important characteristics affected by tillage practices which impact on the physico-chemical 
and hydro-thermal regime in soil, and ultimately crop yield. Additionally, the effect of tillage on the environment 
varies across farms geographically since the impacts of cultivation on soil organic matter and net greenhouse 
balance depends on soil type, climatic variables and management 15 . 

No previous studies have considered the effect of the soil porous architecture created by tillage on net balance of 
greenhouse gas emissions. Traditional methods for inferring soil structure such as soil moisture retention curves 
are limited as they are destructive and do not provide the soil pore size distribution in three dimensions 18 . 
However, imaging technologies such as X-ray Computed Tomography (CT) can be used to reveal the undisturbed 
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Table 1 Selected physico-chemical properties of soils unc 


er zero tillage and conventional tillage 












Depth Bulk density Shear strength Soil moisture 

t;||_,__ t,- m \ IKAi-t /AAP«\ f°/\ 

linage icmj (ivig m } \iv\ro} \/o} 


SOM NH 4 -N (mg 
pn \/o) kg ^11) 


N0 3 -N 
(mg kg _, soi |) 


Microbial C 
(mg kg-'soii) 


Microbial N 
(mg kg" 1 ^1) 


Zerotilled 0-10 1 .16 ± 0.04 25.7 ± 1 .47 31 .29 ± 1 .40 
10-20 ND** ND 27.90 ±1.36 

Tilled 0-10 1.09 ±0.04 12.0 ±1.12 26.98 ± 1 .06 
10-20 ND ND 24.96 ±1.11 


6.98 ±0.1 3 7.81 ±0.44 2.59 ±0.10 
7.32 ±0.10 7.41 ±0.42 2.42 ±0.08 
7.22 ±0.14 6.59 ±0.42 2.51 ±0.16 
7.29 ±0.13 6.15 ±0.40 2.30 ±0.14 


0.66 ± 0.05 
0.45 ± 0.04 
0.62 ± 0.06 
0.54 ± 0.06 


591.8 ± 
442.2 ± 

434.9 ± 
402.5 ± 


55.0 104.9 ± 

26.6 77.3 ± 
44.3 73.4 ± 

39.7 66.6 ± 


7.92 
5.1 1 
5.1 1 
3.79 


*Mean ± Standard Error of mean (n — 33). 
**ND- not determined. 



structure, aggregation and pore characteristics at high resolutions 
(e.g. microscale <100 urn). In this study we sought to evaluate the 
impact of zero tillage and conventional tillage on soil pore character- 
istics, carbon sequestration and GHG emissions. We hypothesised 
that zero tillage improves C sequestration and reduces GHG emis- 
sions compared with conventional tillage through the enhanced 
development of the soil porous network associated with less anthro- 
pogenic disturbance. 

Results 

Soil physical properties. Soil texture varied between the different 
experimental sites ranging from heavy clay soils to lighter sandy soils 
(see supplementary Table 1). Crucially there was no significant 
variation in soil texture between paired fields of conventional and 
zero tilled soils (P > 0.05). Zero tilled soils had a higher bulk density 
(1.16 Mg rrT 3 ) than tilled soils (1.09 Mg reT 3 ) (Table 1, P < 0.001) 
which was not influenced by length of zero tillage management 
which ranged from 5-10 years (P > 0.05). Zero tilled soils had an 
increased shear strength (26 MPa) compared to tilled fields 
(12 MPa) (Table 1, P < 0.001), which was also independent of the 
duration of zero tillage (P > 0.05). Soil moisture content 
(volumetric) was significantly higher under zero tilled soils (29.6%) 
compared to tilled soils (26.0%) (P < 0.01), regardless of the duration 
of zero tillage (P > 0.05). 



Soil pore characteristics. X-ray CT measured soil porosity was 
significantly higher under tilled soil (13.6%) than zero tilled soil 
(9.6%) (P < 0.001, Figure la). The porosity in the surface layer (0- 
10 cm) of tilled soils were 46.9% higher than under zero tilled soils 
and 33.2% higher in tilled compared to zero tilled soils in the 10- 
20 cm layer (P < 0.001). Soil pore size followed similar pattern to soil 
porosity (Figure lb). Pore size significantly varied with tillage type 
and soil depth with increased pore size at the surface layers of tilled 
soil (Table 2, P < 0.05). Pores in tilled soils were twice as large 
(0.52 mm 2 ) as those in zero tilled soils (0.27 mm 2 ) (P < 0.01). The 
largest pore sizes were recorded in the 0-10 cm layer (0.55 mm 2 ) as 
opposed to the 10-20 cm layer (0.24 mm 2 ) (P < 0.001). The surface 
area of the total soil pore system was higher in tilled soils (Figure lc, P 
< 0.001). The surface area of pores was also greater in the 0-10 cm 
depth (1.83 mm 2 ) than the 10-20 cm depth (1.07 mm 2 ) across both 
tilled and zero tilled soil treatments (P < 0.01). 

Soil chemical and biological properties. Zero tilled soils contained 
significantly more soil organic matter (SOM) than tilled soils (P < 
0.001). Soil from the 0-10 cm layer contained more SOM than soils 
from the 10-20 cm layers in both zero tilled (7.8 and 7.4% at 0-10 cm 
and 10-20 cm respectively) and tilled soils (6.6% at 0-10 cm and 
6.2% at 10-20 cm) (Table 1, P < 0.001). There were no significant 
effects for duration of zero tillage on soil organic matter (Table 2). 
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Figure 1 | Soil pore characteristics under zero tilled and tilled managed soil derived from X-ray CT. (a) soil porosity (b), mean soil pore size (c) 
and surface area of soil pores at the surface (0-10 cm) and sub-surface layers (10-20 cm) in zero tilled and tilled soils (average values for different sites and 
standard error of the mean are shown, n = 33). 
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Table 2 Statistical output from linear mixed modelling (texture, tillage, duration, depth) for th 
under zero tillage and conventional tillage (F^fi ,df2) statistic) 


e physicoc 


hemical 


characteristics of soils 


Parameter 


Clay (%) 


Tillage 


Duration of 
non-tillage 


Depth 


Tillage 


x depth 


Duration of non- 
tillage x depth 


Moisture content 


A 07 * 
0.y/(l,58) 


1 7.bO(lJ0) 


ns 


32.2y ( ,,63) 


ns 


ns 


Porosity 


A ~7C\ * 
6./0 (1 ,3 2 )* 


1 0.4V(i ,9) 


ns 


59.3„, 63 )*** 


15.86 


1,63)*** 


ns 


Pore size 


N .J 1 (1,21) 


14.2 1 (i,9|** 


ns 


1 /.Z063 


4.89 


*# 

1,63) 


ns 


Pore area 


14.71(1,36)*** 


1 /.(J 1 (1,9) 


ns 


47.71 ( ,,63)*** 


8.36 


*# 

1,63) 


ns 


pH 


A "70 * 

6.7 46)* 


ns 


ns 


38.49,,, 63) 


15.78 


1,63)*** 


ns 


oUM 


ns 


JJ.z4|i ,9) 


ns 


OA IT *** 

84. 13,,, 63) 


ns 


ns 


NH4-N 


J-°°(l,44) 


ns 


ns 


7-52, ,,63) 


ns 


ns 


NO3-N 


ns 


ns 


ns 


29. 8,, 63)*** 


5.03 


* 

1,63) 


ns 


MBC 


ns 


33.96,, 9,*** 


ns 


37. 14,', ,63)*** 


35.67 


1,63)*** 


A P.O., 


MBN 


ns 


25. 85(, V** 


ns 


20.42,,, 63)*** 


7.44 


1,63) 


ns 


co 2 a 


ns 


8.91(,;, 3 |* 


ns 






co 2 b 


ns 


11.12(1,,,,** 


ns 










CH 4 ° 


ns 


5.79,i,,9)* 


ns 










CH 4 b 


ns 


4.99,i,, 8| * 


ns 










N 2 0° 


ns 


10.04,,,u,** 


ns 










N 2 O b 


ns 


6.38,,,, 4|* 


ns 










Subscripted numbers indicate degrees of freedom for F 

biomass nitrogen, superscripts a and b following CO2 

*"p< 0.001. 

"p<0.01. 

*p<0.05. 


value; dfl = numerator df, df2 = denominator df, ns: non-significant, SOM: soil organic matter, MBC 
CH4 and N2O represents potentials expressed in mg m~ 2 fr 1 and ng g~' h~', respectively). 


microbial b 


omass carbon, MBN: microbial 



Neither ammonium (NH 4 -N) nor nitrate (NO3-N) content in soil 
was affected by tillage. Soil from the upper 10 cm contained signifi- 
cantly higher NH 4 -N than the 10-20 cm layer (Table 1, P < 0.01). 
Nitrate (NO3-N) followed a similar trend to NH 4 -N. Tillage type and 
duration did not influence the NO3-N content (P > 0.05). Soil depth 
significantly influenced NO3-N content (P < 0.001) with highest 
amount in the surface layer (0-10 cm) under both zero tillage and 
conventional tillage. 

Zero tilled soils contained significantly more microbial biomass 
carbon than tilled soils (P < 0.001). The mean microbial biomass 
carbon under zero tilled soil was 517.0 mg kg -1 soil compared with 
418.7 mg kg -1 soil in tilled soils. Microbial biomass carbon was 
significantly higher in the 0-10 cm layer (517 mg kg -1 soil) than 
the 10-20 cm layer (419 mg kg -1 soil) under zero tillage and con- 
ventional tillage (P < 0.001, Table 1). Significantly higher microbial 
biomass carbon was recorded at the 0-10 cm layer in zero tilled soil 
(591.8 mg kg -1 soil) with a significant tillage and depth interaction 
(P < 0.001). However there was no significant effect of duration of 
zero tillage (Table 2). 

Tillage and soil depth significantly influenced soil microbial bio- 
mass nitrogen (Table 1 and 2). Zero tilled soils contained a higher 
microbial biomass nitrogen (91.1 mg kg -1 soil) than tilled soil 
(70.0 mg kg" 1 soil) (P < 0.001). Surface layers (0-10 cm) main- 
tained more microbial biomass nitrogen than sub surface layers 
(10-20 cm) under both zero tilled soils and tilled soils. 

Fluxes of greenhouse gases. Potential C0 2 flux was higher from 
tilled soil than zero tilled soil (P < 0.05, Figure 2a). Potential C0 2 
fluxes under zero tilled soil ranged from 47 to 216 mg m~ 2 h~' with a 
mean value of 141 mg nT 2 h~' whilst under tilled soil it ranged from 
119 to 236 mg m~ 2 h~' with a mean value of 171 mg m~ 2 h -1 . The 
potential C0 2 flux on a per soil weight basis was also higher under 
tilled soil (873 ng g~' h~' soil) compared to zero tilled soil (688 ng 
g" 1 IT 1 soil) (P < 0.01, Figure 2d). 

Potential CH 4 fluxes were generally positive and higher from tilled 
soils (0.044 mg m~ 2 h -1 or 0.22 ng g~' soil) compared to zero tilled 
soil (0.018 mg nT 2 h" 1 or 0.09 ng g" 1 h" 1 soil) (P < 0.05, Figure 2b 
and 2e). In contrast, potential N 2 0 emissions were higher under zero 
tilled soil (0.63 ng g~' h" 1 ) than tilled soil (0.36 ng g~' h" 1 ) (54% 
higher under zero tilled soil when measured on a soil area basis and 



77% on a soil dry weight basis compared to tilled soil) (P < 0.01, 
Figure 2c and 2f). 

The net global warming potential calculated as per IPCC 19 was 
significantly higher from tilled soil than zero tilled soil. Tilled soil 
produced 31% on an area basis or 26% on a weight basis greater 
global warming potential (GWP) than zero tilled soil (P < 0.05, 
Figure 3). There was no evidence to suggest that the different dura- 
tion of zero tillage considered in this study, (5-10 years) affected net 
emissions of greenhouse gases. 

Relationship between greenhouse gas fluxes and soil properties. 

Potenital C0 2 fluxes were predicted by a multiple regression model 
(P < 0.001) including bulk density (BD), microbial biomass carbon 
(MBC) and soil porosity (P) which accounted for 69.9% of the 
variation. The optimal model for the potential C0 2 flux is 
provided in the equation (1). 

C0 2 flux (mgm~ 2 h _1 ) = 124.1 -39. 1BD 

1 ' (1) 

+ 0.0412MBC + 3.689P 

In this model the soil porosity contributed to c. 40% of variation, 
much higher than the individual contribution by any other 
parameter, as illustrated by retaining the parameter when fitting last 
to the model. Together microbial biomass carbon and bulk density 
contributed to 30% of the total variation (Figures 4a, 4b and 4c). 

Only soil shear strength (SS) explained variation (18%) in the 
potential CH 4 flux (Equation 2, Figure 4d, P < 0.01). 

CH 4 flux (mgnr 2 IT 1 ) =0.05344 -0.001078SS (2) 

The optimal model in equation (3) for potential N 2 0 flux accounted 
for 62.0% of the variation and included soil moisture (SM), microbial 
biomass nitrogen (MBN) and microbial biomass carbon (MBC) 
(Figures 4e and 4f, P < 0.001). 

N 2 0 flux (mgnT 2 h~') = - 0.0746 + 0.002057SM- 

(3) 

0.00049 MBN + 0.0003104MBC 

Individually microbial biomass carbon explained the greatest pro- 
portion (20.8%) of the total variation when fitted last in the model. 
Removing soil moisture and microbial biomass nitrogen separately 
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Figure 2 | Fluxes of greenhouse gas from zero tilled and tilled soil, (a) C0 2 expressed in mg C0 2 -C m~ 2 h -1 , (b) CH 4 expressed in mg CH 4 -C m~ 2 h~ 
(c) N 2 0 expressed in mg N 2 0-N rrT 2 h _1 , (d) C0 2 expressed in ng C0 2 -C g -1 h _I , (e) CH 4 expressed in ng CH 4 -C g -1 tr 1 and (f) N 2 0 expressed 
in ng N 2 0-N g -1 h~ (average values for different sites and standard error of the mean are shown, n = 33). 



from the model did not substantially decrease the amount of vari- 
ation explained suggesting that these factors were confounded. 

Discussion 

We have demonstrated tillage practice has the potential to strongly 
influence release of C0 2 , CH 4 and N 2 0, through its impact on soil 
biophysical properties across a wide range of soil textures. However, 
the main driving factors and the direction of change varied among 
the three GHGs measured. The higher C0 2 release found in response 
to tillage highlights the role of ploughing in the breakdown of soil 
aggregates and exposure of organic materials for microbial decom- 
position 20 . Soil pore characteristics (previously ignored in similar 



studies), such as total porosity and pore size, were a stronger pre- 
dictor of C0 2 flux than soil organic matter and microbial biomass 
carbon, which has not been previously reported. The effects of zero 
tillage was to reduce soil porosity by 33%, which lead to 21% reduc- 
tion in potential C0 2 efflux. These results demonstrate the increased 
soil porosity under conventional tillage favours the respiration of 
aerobic organisms by improving movement of water and air through 
the soils 21 with important implications for C0 2 emissions. In parallel, 
strong effects of soil bulk density on C0 2 production from soil cores 
have been shown by Beare et al. 22 who found 2.3 times more C0 2 
production under uncompacted soil than in compacted soil. The 
potential C0 2 flux data presented here (47 to 235 mg irT 2 h -1 ) is 
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Figure 3 | Global warming potential under zero tilled and tilled soils. (Average values for different sites and standard error of the mean are shown, 
n = 33). (a) GWP expressed in terms of mg rrT 2 tr'and (b) GWP expressed in terms of ng g _1 h~'. 



SCIENTIFIC REPORTS | 4:4586 | DOI: 10.1038/srep04586 



4 



a) 



250 



200 



| 150 
« 100 



•Zero tilled OTilled 



O 50 y=-167.93x + 344.87 # « 
U R J = 0.40 



0.6 0.8 1 1.2 1.4 

Bulk density' (Mg nv 5 ) 



1.6 



b) 

250 

% 200 

| ,50 

§ 100 
C 

8 50 



•Zero tilled OTiUed 



t 



y =0.134x+ 116.78 
R ; = 0.30 



500 1000 1500 

Microbial biomass carbon (mg kg 1 ) 





1.0314 



Soil moisture content (%) Microbial biomass carbon (mg kg 1 ) 



Figure 4 | Illustration of important relationships between soil biophysical properties and GHG release, (a) soil bulk density and C0 2 flux from soil; F lj64 
= 42.08, P < 0.001 (b) microbial biomass carbon and C0 2 flux; F 1>64 = 5.89, P < 0.05 (c) soil porosity and C0 2 flux; F 1>64 = 110.14, P < 0.001 
(d) soil shear strength and CH 4 flux; F lj64 = 14.08, P< 0.001 (e) soil moisture content andN 2 0 flux; F 1>64 = 12.62, P< 0.001 and (f) microbial biomass 
carbon and N 2 0 flux; F 1>64 = 69.5, P < 0.001. 



in the range of that reported from laboratory incubations of soils 
from 13 European sites including, arable land (47 mg mT 2 h -1 ) and 
grassland (186 mg m~ 2 h -1 ) 23 . Similar effects of tillage on C0 2 fluxes 
were found by Ball et a/. 24 investigating in situ C0 2 fluxes, they 
attributed the greater C0 2 efflux to the larger pores created by tillage. 
Potential CH 4 flux ranged from 0.0025 to 0.16 mg m~ 2 h^ 1 , which is 
high compared to values reported by Schaufler et al. 23 : e.g. average 
CH 4 flux in arable land was 0.0014 mg m~ 2 h -1 and in grassland it 
was 0.0005 mg m~ 2 h -1 . Despite the less porous and wetter status of 
zero tilled soils, which normally promote CH 4 production 31 , the 
opposite was the case here which may be due to increased activity 
of methanotrophic bacteria 32 . The reduced potential CH 4 flux under 
zero tillage was best predicted by soil shear strength which reflects 
the reduced porosity and high bulk density in zero tilled 
soils 17,25,26 .Increased bulk density in soil can prevent flow of CH 4 in 
soil and the resulting enhanced retention of CH 4 in soil may improve 



oxidation by methanotrophs 27 resulting in lower CH 4 emissions. 
Furthermore, the development of methanotrophic populations is 
negatively affected by tillage 28 which are slow to recover 29,30 . The 
potential N 2 0 fluxes measured were comparable to field measure- 
ment by Regina et al. 33 in Finnish soils after 5-7 years of zero till 
management (0.003 to 0.23 mg m" 2 h" 1 ) with significantly higher 
N 2 0 fluxes under zero tillage. They reported 21 to 86% higher N 2 0 
flux in zero till soils when compared to tilled soils. The average 
increased emission of in situ N 2 0 flux under zero tilled soils obtained 
by Oorts et al. 34 was 39% for a 30 year experiment. As with CH 4 , N 2 0 
is produced under reducing conditions in waterlogged and poorly 
aerated soils 35-36 , so we attribute the increased potential N 2 0 emis- 
sions from zero tilled soils in part to the wetter and denser soils found 
under this management regime. In contrast to the potential C0 2 and 
CH 4 fluxes, the potential production of N 2 0 was most strongly 
related to the soil microbial biomass. The greater total soil microbial 
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biomass found under zero tillage may hence play a very important 
role in N 2 0 release. One important aspect of zero tillage is enhanced 
crop residue retention resulting in greater SOM content. Given the 
importance of an adequate supply of labile substrates for the denit- 
rifying bacteria 35 , it may also be that the crop retention under zero 
tillage drives greater N 2 0 release. 

Considering the GHGs together, tilled soil produced 20% greater 
net global warming than zero tilled soil indicating a potential for zero 
tillage system to mitigate climate change after only 5 to 10 years since 
conversion (earlier than this was not measured here). In parallel with 
this Del Grosso et al. 37 also reported a 33% reduction in global warm- 
ing potential under zero tillage (0.29 Mg C ha" 1 y" 1 ) compared with 
tilled soil (0.43 Mg C ha" 1 y" 1 ) for major non-rice cropping systems 
in US based on simulation using DAYCENT ecosystem model. Also 
in subtropical conditions, zero tillage has been found to reduce GWP 
by c. 20% 38 . 

Zero tilled soils had enhanced SOM, microbial biomass carbon 
and nitrogen. Importantly, the time during which the soils had been 
under conservation tillage did not influence the SOM content in the 
soil (although only changes between 5 and 10 years were measured), 
suggesting that increases in SOM occurred within five years follow- 
ing conversion to zero tillage. However West and Post 5 in similar 
work recorded a large increase in soil between 5-10 years. The time 
required to reach a steady state for carbon sequestration will vary 
with respect to climate, soil types and the management practices 
followed 39 . 

A very important question that remains to be addressed is how the 
impact of the change in the soil porous architecture brought by 
tillage/zero tillage on net GHG release and the GWP varies spatio- 
temporaly across a greater range of soils types, crops and climate than 
those explored in our study. With reduced tillage practices becoming 
more prevelant globally, it is important to further understand the 
impacts of this on the biophysical evolution of the soil environment 
at both micro and macroscales. It is clear from this study that the 
modification of soil structure by tillage plays a crucial role for GHG 
release. Our study was based on analysis of intact cores removed 
from the field. To fully account for the impact of zero tillage on 
GHG release it is important to extend this work to in situ field 
measurement through the year to account for variation in weather 
and crop development. In conclusion, we have shown soils under 
zero tillage increased potential N 2 0 emissions, but this is counter- 
balanced by a significant reduction in potential C0 2 and CH 4 emis- 
sions which is closely linked to the geometry of the soil porous 
architecture. To evaluate the potential of zero tillage as a tool for 
mitigation of climate change, there is a need to further assess its 
impact on yield to ensure a balance between climate change mitiga- 
tion and food security is achieved. 

Methods 

Soils. A selection of 22 farms from Leicestershire, Nottinghamshire and Lincolnshire 
in the East Midlands of the U.K. were chosen for this study. All sampling sites 
comprised pairs of intensely tilled farms and farms where zero tillage practices were 
followed, located directly adjacent to each other. The zero tilled soils had been 
managed in this way for a minimum of 5 years to a maximum of 10 years. In fields 
under zero tillage, stubble was left at the surface after harvest of the previous crop. 
Seed drilling was carried out between the root stocks of previous crop using a range of 
min-till seed drills. The crops cultivated under zero tilled and tilled sites were wheat, 
oil seed rape and oats. The tilled soil sites were annually ploughed to depths of 
20-25 cm and contained the same crops as the zero tilled fields at the time of 
sampling. Sampling was undertaken shortly after seedbed preparation and sowing so 
as to minimize any effect of the emerging root system on soil structure. 

Intact soil cores were collected using a manual core sampler that used transparent 
sample liner tubes (Van Walt Ltd, Haslemere, UK). The core sampling was performed 
to a depth of 20 cm with a diameter of 5 cm and in triplicate. The samples were 
labelled and sealed in plastic bags before transporting to the laboratory. Samples were 
stored at 4°C until measurements were taken (<2 weeks). Bulk soil samples of about 1 
kilogram were also collected from two depth ranges (0 to 10 cm and 10 to 20 cm) and 
stored at 4 C until measurement. Smaller soil cores were collected in the field using 
stainless steel cylinders (radius 3.4 cm, height 4 cm) for measurement of bulk 
density 40 . 



Soil physical properties. Soil shear strength was recorded in the field using a Pilcon 
120 kPa hand vane from the upper 50 mm of soil. Similarly the volumetric water 
content of the surface layer of soil (0-10 cm depth) was recorded using a Delta-T 
Theta probe connected to a Theta meter. All observations were recorded in triplicate 
for each field. Particle size analysis was performed using the hydrometer method 41 . 
Soil textural classification was made according to European classification using 
60 um as the upper limit for silt 42 . 

X-ray Computed Tomography (CT). Prior to the study of GHGs, the soil core 
samples were subjected to morphological analysis using an X-ray CT scanner 
(Nanotom, Phoenix X-ray, GE Sensing and Inspection Technologies GmbH, 
Germany) to visualise and measure the internal soil structure. The cores were scanned 
at a voltage of 140 kV and a current of 100 mA. A copper filter of thickness 0.25 mm 
was used to minimise artefacts such as beam hardening. The image resolution was 
64 um per voxel. The soil core was positioned vertically onto the scanner platform. 
Each scan lasted 100 minutes per core, scanning both top and bottom 10 cm portions 
in a split scan. Whilst it is possible to achieve much faster scan times than this, a larger 
scan time was used to achieve the highest possible image quality. For each scan 1000 
images were collected. The obtained images were visualised using the software, VG 
StudioMax (Volume Graphics). The images were converted to the tiff format and 
analysed using ImageJ 43 to study the soil pore characteristics. A rectangular region of 
interest (27.94 X 27.94 mm 2 ) was selected to avoid the edges of the soil cores. In 
addition the first 100 images each from the beginning and end of the scan were 
discarded due to cone beam artefacts. The images were sharpened to highlight the 
image features and then smoothed by a median filter before being converted to the 
binary scale using the minimum threshold algorithm in ImageJ. Both dark and bright 
outliers were removed and the 'fill holes' function was used to minimise noise. 
Measurements on soil physical features were obtained on the binary images which 
included porosity, number of pores, pore size and surface area of pores (Figure 5). 

Soil chemical and biological properties. Soil pH was determined on air dried 2 mm 
sieved soils using 1 : 2 soil to water ratio using a combined glass electrode. Total soil 
organic matter (SOM) content in soil was determined by loss on ignition following 
igniting oven dried soil at 550°C in a muffle furnace. For the measurement of 
ammonium and nitrate (NH 4 -N and N0 3 -N) concentration, 6 g of field moist soil 
was used. An extraction was carried out using 40 ml of 2 M KC1 by shaking and 
filtration. Ammonium in the extracts was determined colourimetrically 44 . A suitable 
aliquot of the filtrate ( 1 ml) was made to react with phenol and hypochlorite to form a 
blue indophenol complex in solution. The concentration of ammonium in solution 
was measured by comparing the absorbance with known standards prepared using 
NH 4 C1 at a wavelength of 635 nm. For the determination of N0 3 -N, nitrate in a 
suitable aliquot of KC1 extract was reduced to nitrite using spongy cadmium, which 
was further complexed to form a red azo-species in solution. The concentration of 
NO3-N was measured by comparing the absorbance with known standards of KN0 3 
at a wavelength of 543 nm 45 . Field moist soil samples (both surface 0-10 and 
subsurface 10-20 cm depths) were used for the estimation of microbial biomass 
carbon and nitrogen by the chloroform fumigation- extraction technique 46 . Samples 
were incubated in the chloroform environment in presence of soda lime. The 
extraction was carried out using 0.5 M K 2 S0 4 at the start of fumigation in un- 
fumigated samples and 24 hour after fumigation in fumigated samples. Microbial 
biomass carbon and nitrogen in the extracts were analysed using a Shimadzu CN 
analyser (TOC-V CPH Shimadzu). The results were corrected using the value of 0.45 
for both carbon and nitrogen 47 . 

Potential fluxes of greenhouse gases. Cores were removed from the 4°C 
environment and kept at a constant temperature of 16°C for 48 hours to activate and 
stabilise the biological activity. Gas sampling was performed by placing cores in 1.5 
litre plastic jars (20 cm height and 10 cm diameter) with a septum on the top to aid 
gas sampling using a 20 ml syringe. The air in the headspaces was mixed, before 
sampling at time intervals 0, 15, 30 and 60 minutes using 20 ml syringes. The 
collected gas samples were stored in airtight pre-evacuated glass vials and analysed for 
concentration of C0 2 , CH 4 and N 2 0 using gas chromatography equipped with a 
Thermal Conductivity Detector (TCD), Flame Ionization Detector (FID) and an 
electron capture detector (ECD) (GC-2014, Shimadzu). The fluxes of these samples 
were calculated using linear regression of the gas concentration against time. The 
GHG data was converted to mass per volume and mass per weight basis by the use of 
ideal gas equation and the molecular mass of each gas 48 . 



Where n is the number of moles of C0 2 , N 2 0 or CH 4 , P is atmospheric pressure 
(~1 atm), V is the volume of head space (dm -3 ), R is the ideal gas constant 
(0.08205746 L atm K" 1 mol" 1 ) and T is the temperature of sampling (273.15 + room 
temperature in °C). From this the flux of gas was measured. 

nm 

E= — X 1000 (5 
at v 1 

Where E = flux of each gas in mg m" 2 h~ ! , n — number of moles of C0 2 , N 2 0 or CH 4 , 
m = molar weight of C0 2 (44.01), N 2 0 (44.01) or CH 4 (16.04), a = area of the soil 
core used and t is the time in hour. Finally total greenhouse balance or net global 
warming potential (GWP) was calculated in C0 2 -equivalents 19 using the following 
equation. 
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Figure 5 | Non-destructive 3-D imaging of soil by X-ray CT. Examples of Tilled (A-C) and Zero Tilled (D-F) soils. (A&D): 3D rendered grayscale 
density map of soil cores showing a virtual 'cut-out' to the revealing clear differences soil structure between the two soils. (B&E): Thresholded 3D image 
highlighting 'solid' soil in brown and 'pore' space in white. (C&F): Visualisation of pore space only highlighting high connectivity of pores in the tilled 
soils and the presence of numerous biopores in the zero tilled soil. Scale bar =10 mm. 



GWP = 



C0 2 x 44 
12~ 



+ (CH 4 x23) + (W 2 Ox296) 



(6) 



Statistical analysis. Each site consisted of a pair of fields; one of which was ploughed 
and the other had been tilled for a number of years. The sites were in areas consisting 
of a range of soil types in different geographical regions although at each site the tilled 
and zero tilled fields were located adjacent to each other. Samples were taken at a 
number of random locations in each field and at two soil depths (0-10 and 
10-20 cm). The variation in soil properties in response to tillage and soil depth was 
analysed as split-split plot design in a linear mixed model with site, field and location 
within fields as random effects. Tilled vs non-tilled, soil depth and their interaction 



were considered as fixed effects. The variation among just the zero-tilled fields was 
further partioned to test for a trend in response to the number of years since adoption 
of zero tillage. This test was thus orthogonal to the tilled vs zero-tilled contrast as was 
its interactions with soil depth. To account for potential differences with respect to 
soil texture, the clay content of the soil was considered as a covariate by including it as 
a fixed effect in the model. The covariate could account for significant amounts of the 
random variation among fields and locations. In such cases, by reducing the 
unexplained residual variation, the model including the covariate is likely to be more 
sensitive for detecting tillage effects than a model without the covariate. Multiple 
linear regressions were used to predict the best model describing the fluxes of GHGs 
from soil. The maximal model consisted of all the physical, chemical and biological 
properties studied in this experiment. By using a stepwise backwards elimination 
process, only the variables that contributed significantly to the model and reduced the 
residual sum of squares were retained in the model. For illustrative purposes we also 
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carried out the single linear regression between the parameters that contributed to the 
multiple regression models. All tests were performed using Genstat {14 th Edition, 
VSN International Ltd, UK). 
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